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SUMMARY 

A program was designed for the calculation of the performance, weight, and 
area characteristics of a single-panel central-fin-tube flat-plate sensible-heat 
radiator for a set of thermodynamic and fluid mechanic conditions. The 
FORTRAN 7094 main program and its subroutines are reported herein. The program 
is based on the analysis, equations, and techniques developed and reported in 
NASA Technical Note D-2839. 

A typical computer printout sheet is included and a discussion of the def- 
initions of the FORTRAN symbols and groups of numbers that designate the inputs, 
the error tests, and the outputs is presented. 

INTRODUCTION 

A design analysis and general characteristics of flat-plate central-fin- 
tube sensible -heat space radiators were presented in reference 1. It was shown 
therein that a combination of equations based on the relations of heat transfer, 
fluid mechanics, and meteoroid protection resulted in specific relations for 
calculating radiator geometry, weight, and panel planform area. In order to 
perform the parametric studies reported in reference 1, these equations and re- 
quired inputs were programed for an IBM 7094 electronic digital computer. 

The purpose of this report is to describe the FORTRAN 7094 main program 
and its subroutines. Some phases of the computational procedure (e.g. , 

Kalaba's method) will be discussed in more detail. Computer program diagrams 
are included to illustrate the procedure and to show how the various elements 
and subroutines are interconnected to generate the required outputs. 

In addition, a typical printout sheet will be presented and explained. It 
contains the inputs and the outputs for two sensible-heat radiator examples: a 

Brayton cycle case and a secondary cooling radiator case. Various FORTRAN 
symbols and groups of quantities without symbols will be defined as they are 
written and appear on the printout sheet. The definitions correspond to the 
definitions of the symbols used in reference 1. 




SYMBOLS 


A area, sq ft 

D diameter, ft 

F angle factor 

h heat -transfer coefficient, Btu/(sec)(sq ft)(°R) 

L half -fin width, ft 

Re Reynolds number 

R 0 outside tube radius, ft 

T absolute temperature, °R 

t fin thickness, ft 

U G overall heat-transfer coefficient, Btu/(hr)(sq ft)(°R) 

W weight, lb 

Z tube length, ft 

X. heat -transfer parameter 

& armor thickness, ft 

t] fin-tube effectiveness 

9 dimensionless temperature ratio 

\ dimensionless conductance parameter 

Subscripts : 

H header 

j at location j 

o outside 

s sink 

t tube 

tot total 

X at location X 

COMPUTATIONAL PROCEDURE AND COMPUTER PROGRAM LISTING 

Procedure 

As shown in the analysis section of reference 1, a combination of equations 
based on the relations of the heat transfer, fluid mechanics, and meteoroid pro- 
tection resulted in certain relations for calculating radiator size, weight, and 
panel planform area. In order to perform the parametric studies, these equations 
and required inputs were programed for an IBM 7094 electronic digital computer 
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The FORTRAN 7094 program consists of two versions of the main program and 
four subroutines. The only difference in these versions is that the first ver- 
sion finds the fin-tube effectiveness by solving the proper equations using 
subroutines TOJ, DEQ2, and FXX. The second version using subroutine TABLE finds 
fin-tube effectiveness by interpolating in a table listing fin-tube effective- 
ness , which is read into the machine. Hereinafter, the main program refers to 
either version. 

The flow diagram of the main program is shown in figure l(a) . The main 
program is written for either a gas, a liquid, or a liquid metal as the working 
fluid. For the radiator using gas as the working fluid, the program involves 
the ratio of header surface area to tube surface area A^/A^ . For liquids and 
liquid metals the ratio A H /At does not appear, since the header area is 
assumed to be negligible compared with the tube surface area. 

Block 1 of figure l(a) contains the main program inputs. They are the 
physical and thermal properties of the working fluid and the radiator structural 
materials (such as, sonic velocity, viscosity, specific heat, density, thermal 
conductivity, emissivity), as well as meteoroid protection parameters (meteoroid 
flux, density, operation time, probability of no meteoroid penetration, occlu- 
sion and spalling factors). Also included are the flow rate, inlet and outlet 
temperatures, the pressures and pressure drops, the accuracy percentages of the 
computations, the branching command numbers, and the number of points in the 
various meshes. 

Parts of some equations can be computed based on the inputs and constants 
listed in the previous paragraph. This is done in block 2 of figure l(a) . The 
next step is to read and write the input parameters: inside tube diameter, fin- 

tube profile ratio, and fin conductance parameter at the radiator inlet. There 
are branches throughout the program that select either gas, liquid, or liquid 
metal radiator equations, variables, and constants. The program is written for 
tubes and headers with inside liners. Both turbulent -flow and laminar-flow 
equations are contained in the program for all three mediums: gas, liquid, and 

liquid metal. First, a solution is obtained for given inputs from turbulent- 
flow equations. After the solution is converged, the Reynolds number is check- 
ed. Then the Reynolds number is used as a test either to read the next set of 
inputs (if the Reynolds number indicates turbulent flow) or to find a second 
solution from the laminar-flow equations. A single result is printed out when 
the Reynolds number is less than 2300 or greater than 3000. If the Reynolds 
number is between 2300 and 3000, both laminar and turbulent flow results are 
printed. 

The noteworthy feature of the main program is the iteration for the over- 
all heat-transfer coefficient. The iteration starts with a guessed value and 
proceeds by using the method of false position (Regula Falsi) that is based on 
the principle that a curve over a short interval can be approximated by a 
straight line between two values, each being on opposite sides of the right 
value. The same method is used to calculate the ratio Ag/A^. for the gas work- 
ing fluid. The overall heat-transfer coefficient iteration is accomplished 
within the Ag/Aj. iteration. 

During both iteration processes the values for the tube and header wall 
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(a) Main program. 

Figure 1. - Computer program flow diagrams.. 
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(b) Subroutine TOJ. 
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(c) Subroutine FXX. 

Figure L - Continued. Computer 
program flow diagrams. 


thickness , the tube outside diameter , the tube length, 
number of tubes^ tube and liquid content weights, and 
the panel planform area are also obtained. With this in- 
formation, the half -fin width, the fin thickness, fin 
weight, and header dimensions and weights are obtained. 


Within the iteration of overall heat-transfer coef- 
ficient, the program solves for the tube surface tem- 
perature of each isothermal elemental strip by Newton’s 
method. The first version requires a subroutine TOJ, 
the flow diagram for which is shown in figure l(b). The 
fin-tube effectiveness is also computed by the subroutine 
TOJ. This involves solving the second-order differential 
equation 



(la) 


where 


F X “ F X-1 + F X-2 


0=1 at X = 0 

ae/ax = o at x = l 


(it) 


( 2 ) 


ty using subroutines DEQ2 ana FXX. These aaaitional sub- 
routines will be aiscussea later. As part of Newton's 
metho!, the lerivative aT) is compute!, ana a value of 
T 0j j is obtaine!, which is checked against the guessed 
value aepictea in the first block (fig. 1(b)). The iter- 
ation is repeatea until the new value of T 0 ,j agrees 
with the previous value within the require! accuracy. 

The subroutine FXX (fig. 1(c)) computes the angle 
factors between fin an! tube at each mesh point for use 
in the subroutine DEQ2 (fig. 1(a)). Subroutine FXX also 
obtains integrals in!epen!ent of 0 to be use! in sub- 
routine TOJ (fig. 1(c)) for calculating the integral in 
the expression for the fin-tube effectiveness (ref. l) . 
The subroutine DEQ2 is use! to solve the secon!-or!er 
aifferential equation (eqs. l(a) an! (b)) with its 
boun!ary conaitions (eq. (2)). The method use! in this 
program was devised by Kalaba (ref. 2). With this 
metho!, 0 4 is approximate! by a linear function 0, 
namely, by the first two terms of its Taylor series. The 
central aifferences are applie! to a^/aX 2 . The linear 
approximation, the aifference approximation to ! 2 e/dX 2 , 
and the boun!ary conaitions are use! to obtain a set of 
linear equations whose correspon!ing tri-!iagonal matrix 
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Start 
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(d) Subroutine DEQ2. (e) Subroutine TABLE. 

Figure 1. - Concluded. Computer program flow diagrams. 


is reduced to two-diagonal form. The temper- 
ature Q profile is then solved hy backward 
substitution. An initial 6 profile guess 
is used and each successive iteration yields 
approximately one significant figure of 
accuracy. 

By using the second version of the main 
program, computation time can be signifi- 
cantly reduced with the introduction of a 
table. The subroutine called TABLE 
(fig. 1(e)) now is used instead of sub- 
routines TOJ, FXX, and DEQ2. 

The values of the fin-tube effective- 
ness were first obtained from subroutines 
TOJ, FXX, and DEQ2 by using the fin-tube 
profile ratio, conductance parameter, and 
sink temperature ratio as independent vari- 
ables. These values are then tabulated for 
use in subroutine TABLE. The first version 
may be used if these parameters are out- 
side the range of TABLE or if it is desired 
to avoid the slight inaccuracy caused by 
using linear interpolation of nonlinear 
functions. 
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Program Listings 


$ I BFTC GASLMC DECK 

C SAULE - GAS » L I OU ID*LIOUID METAL RADIATOR 0010 

DIMENSION T ABL (13*36)* TWR ( 1 00 ) » DEL ZNW ( 1 00 ) * DELZ W ( 1 00 ) » 0020 

1 FLWRIT(IOC) , 

2 DI J( 20) ,FLRJ (20) ,FLAMRJ ( 20) 

COMMON TABL.FLR.ENC.ETATOT *SLPME*TUSE 

READ ( 5, 130) ( (TABL ( I ,J) , J=1 ,36 ) , 1 = 1,13) 0060 

130 FORMAT ( 12F6. 3 ) 0070 

132 FORMAT ( 9F6. 2 ) 0090 

READ < 5.134)PIE»FNMESH,FTEST,TEST,TESUH*FATST 0100 

134 FORMAT ( El 3.6.6E10.3 > 0110 

WRITE (6,467) 0120 

215 READ ( 5.3 ) STEPH, ALPHA, BETA, EMACH.RHAP ,VAVG,PNAN» FCCLU.TDAY, 0130 

1EPSL»ATICK*GSRLC*QLMR*TZWR 0140 

NMESH=FNMESH 0150 

3 F0RMAT(6E12.5) 0160 

4 FORMAT ( PEI 5,5) 0190 

0200 

BRANCH - GAS OR L I QU I D , L I OUI D METAL RADIATOR 0210 

IF GAS CASE, THEN GSRLQ = 0. AND QLMR = 

IF LIQUID CASE, THEN GSRLQ = 1 AND QLMR = 0 
IF LIQUID METAL CASE , THEN GSRLQ = 1 - AND QLMR = 1 


0220 

2 IF(GSRLQ) 160,160,150 

150 READ (5»9)SCVEL*FMAS,VISC»UHETRD»CPH»TINLT,TEXIT,RhAIN,RhAF, 0240 

1 RHAH,RHAT,PDRPT,PDRPH*THERKF,TCGAS,TS»RHAC 0250 

IF (QLMR) 664,664,665 

664 WRITE (6,662) 

GD TO 667 

665 WRITE (6,666) 

667 WRITE (6,4)SCVEL»FMAS,VISC,UHETRD,CPH,TINLT,TEXIT, RHAIN.RHAF, 0260 

1RHAH,RHAT,PDRPT,PDRPH,THERKF,TCGAS,TS,PHAC,STE»H, ALPHA, BET A, 0270 

2EMACH,RHAP,VAVG,PNAN,FCCLU,TDAY,EPSL,ATICK.,FTEST, TEST.TESUH, 0280 

3FNMESH, GSRLQ, QLMR, TZWR 0290 

419 TCLQ=T CGAS 0350 

RHAL=RHA I N 0360 

155 FF= 1 . 0 0370 

GD TO 165 0380 

160 READ (5*9)SCVEL»FMAS,VISC*UHETRD*CPH,TINLT,TEXIT, RGAS.PINLT, 0390 

1 RHAF ,RHAH,RHAT»PTUBT,PTUBH,THERKF,TCGAS»TS,RHAC 
WRITE (6,663) 

WRITE (6,4)SCVEL,FMAS,VISC,UHETRD,CPH,TINLT,TEXIT, RGAS.PINLT, 0410 

lRHAF.RHAH.RHAT.PTUBT.PTUBH.THERKF.TCGAS.TS.RH AC. STEPH. ALPHA. 

2 BET A, EMACH ,RHAP »VAVG.PNAN*FCCLU.TDAY»EPSL*ATICK»FTEST .TEST « 

3 TESUH.FNMESH.FATST, GSRLQ, QLMR, TZWR 

PDRPT=PINLT*PTUBT 0450 

PDRPH=PINLT*PTUBH 0460 

PEXIT=PINLT-( PDRPT+2.0*PDRPH) 0470 

RHAIN=PINLT/(RGAS«TINLT) 0480 

163 RHAUT=PEX I T/(RGAS*TEXIT) 0490 

T I NTD= 1 • 0/T I NLT 0500 

TEXTD=1.0/TEXIT 0510 
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TINTD2=TINTD**2 

0520 


TEXTD2=TEXTD**2 

0530 


FF2=TEXTD2-TINTD2 

0540 


FF3=TEXTD2*TEXTD-TINTD2*TINTD 

0550 


FF=1.5*FF2/(FF3*TINLT) 

0560 

165 

FLTN= 2.0 



DELT = ( T I NLT -TEX I T ) /FNMESH 
TEMPRD=TI NLT-DELT /2 .0 

0580 


HI1=RHAIN/VISC 

0590 


HI 3 = VI SC*CPH/TCGAS 

0600 


IF(OLMR) 204,20a, 201 

0610 

201 

HI4=HI 3**0. 4 

0620 


HI 5 = 0 • 62 5*TCGAS*H 14 

0630 


GO TO 206 

0640 

204 

H I 4=HI 3**0. 3 

0650 


HI 5 = 0. 023*TCGAS*H 1 4 

0660 

206 

DELZN1=EMAS*CPH*DELT 

0670 


DELZN2=3600.0*DELZN1/PIE 

0680 


FMAS2=FMAS*FMAS 

0690 


DELT1=(62.45*RHAP/RHAT)**0.5 

0700 


DELT2=(V AVG/ SC VEL)** (2.0/3. 0) 

0710 


DELT3= ( 0.6747E-04/RHAP)**( 1. 0/3.0) 

0720 


BETAD=1.0/(3. 0*BET A ) 

0730 


DELT4=-TDAY*ALPHA/AL0G(PNAN) 

0740 


DELT5=2.0/(3.C*EMACH*BETA+2.C ) 

0750 


DELT6=PIE*DELT4*DELT5 

0760 


DELTA1=DELT1*DELT2*DELT3 

0770 

5 

DELTA2=2.0*ATICK*FCCLU*DFLTA1 

0780 


FMAS23=20.3/FMAS 

0790 


ZBG3=RHAIN*PDRPT/FF 

0800 


ZBG4=ZBG3**( 5.0/14.0) 

0810 


ZBG5=( 1.0/VISC)**( 1.0/14.0) 

0820 

6 

ZBG7=ZBG4*ZBG5 

0830 


VNE1=4.0*FMAS/(PIE*RHAIN) 

0840 


WT1 =P I E*RHAT 

0850 


TS4 = T S**4 

0860 

7 

WF1 =2. 0*RHAF 

0870 

8 

READ (5,524) JD I , ( D I J ( J ) * J= 1 , JD I ) 

READ( 5,524) JFLR , ( FLR J ( J ) » J = 1 , JFLR ) 

READt 5,524) JBRU* ( FLAMR J ( J ) , J = 1 ,JBRU) .READ 
DO 525 JB 1=1, JBRU 
FlAMRD= FLAMRJ(JBI) 

DO 526 JB2=1 , JFLR 

FLR= FLR J ( JB2 ) 

DO 120 JB3= 1 , JD I 
DI = DI J ( JB3 ) 

WRITE (6,10) DI, FLR , FLAMRD 


10 

FORMAT ( 1 HC , 5H DI=E13.6,3X, 

15H FLR=E13.6,3X ,7HFL AMRD=E1 3 . 6 ) 

0910 



I F ( RHAC ) 460 ,460,461 


460 

DLCC=0.0 



GO TO 464 


461 

D 104= ,04*DI 



IF(DI04-. 00125) 417,415,415 

0310 

415 

DLCC=DI 04 

0320 


GO TO 464 


417 

DLCC=. 00125 

0340 
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464 DLCC2= DLCC+DLCC 
JL23 = 0 
LAOTR= 2 
DEBUG DLCC »D I 


605 IF(GSRLQ) 306.306.310 

306 JAT=1 0940 

JAT 1=0 0950 

JAT2=0 0960 

JAT3=2 0970 

JAT 4 = 1 0980 

JAT5=0 0990 

JAT 6=0 1000 

AHT=0.3 1010 

AHT1=0.3 1020 

308 AHT2=0.3 1030 

GO TO 312 10A0 

310 JAT = 0 1050 

312 GAM 1=STEPH*FLTN*EPSL/P I E 1060 

TLITT1=STEPH*FLTN*EPSL/THERKF 1070 

FNC=FLAMRD 1080 

CALL TABLE 1090 

ETARD=ET ATOT 1100 

DI2=DI*DI 1110 

ZBG6=DI**(12. 0/7.0) 1120 

GAM2=ETARD*I 1.0+FLR) 1130 

11 GAM3=GAM2*GAM1 11A0 

UHETP=UHETRD 1150 

<CNT = 0 H60 

14 UHET =UHETP 11^0 

TEMP1=TEMPRD 1180 

FLAM1=FLAMRD 1190 

ETA1=ETARD 1200 

GAMMA1*GAM3/UHET 1210 

GAM5AV=GAMMA1 1220 

<PP=0 1230 

TDMY=TEMP 1 12A0 

17 TDMY2=TDMY*TDMY 1250 

TDMY4=TDMY2*TDMY2 1260 

FNWT1=GAMMA1*(TDMY4-TS4) 1270 

18 FNWT = FNWT 1+TDMY-TEMP1 1280 

I F ( ABS(FNWT)-FTFST) 25,25.20 1290 

20 IFIKPP-20) 23.23,21 1300 

21 WRITE (6.22) 1 3 1° 

22 FDRMATt 1H0.32H TROUBLE-LOOK AT STATEMENT NO. 20) 1320 

GO TO 120 1330 

23 KPP=KPP+1 1 3 AO 

TDMY3=TDM\2*TDMY 1350 

DNWT =4.0*GAMMA1*TDMY3+1 • 1360 

24 TDMY=TDMY-FNWT /DNWT 1370 

GO TO 17 1380 

25 T SAVE = TDMY 1390 

T2ERC=TDMY**3 l 400 

DELZN=DELZN2/ (TEMP 1-TSAVE) /UHET 1A10 

DZN=DELZN i42 0 

TWR ( 1 ) =TDMY 1^30 

DELZNW ( 1 ) =DELZN 1440 

JEND=NMESH 1450 
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on nn 


FLWRITI1 )= FLAMRD 

30 DO 55 J = 2 * JEND 1 4 60 

KNN = 0 1 47 0 

TEMP2=TEMP1-DELT l 489 

TSAVE3=TSAVE**3 1 49 ° 

DFNTW1=3.0*GAMMA1*FLAM1 1500 

DFNTW2=DFNTW1/(ETA1*TSAVF3) 1510 

TDMY=TEMP2 1520 

35 FLAM2=FLAM1*(TDMY/TSAVE)**3 1530 

36 FNC=FLAM 2 l 549 

CALL TABLE 1550 

37 ETA2=ETAT0T 1569 

GAMMA2=GAMMA1*(ETA2/ETA1 ) 1570 

TDMY2=TDMY*TDMY 1580 

TDMY4=TDMY2*TDWY2 1590 

FNTW1=GAMMA2*(TDMY4-TS4) 1500 

39 FNTW=TEMP2-TDMY-FNTW1 1510 

I F( ABS< FNTW) -TFST )45 ,45 *40 1520 

40 I F t KNN-20 ) 43,43,41 1630 

41 WRITE (6,42) 16A0 

42 FORMAT ( 1H0,32H TROUBLE-LOOK AT STATEMENT NO. 40) 1650 

GO TO 120 1660 

43 KNN=KNN+1 167 0 

TDMY3=TDMY2*TDMY 1680 

DFNTW3=(TDMY4-TS4)*TDMY2 1690 

DFNTWT=DFNTW2*DFNTW3*SLPME 1 7 00 

DFNTW4=4.0*GAMMA2*TDMY3 1110 

DFNTW=-( 1 • 0+DFNTW4+DFNTWT ) 1 ? 20 

44 TDMY=TDMY-FNTW/DFNTW 1 ?3 0 

GO TO 35 l 740 

45 TSAVE=TDMY 1250 

TEMP1=TEMP2 17 60 

GAMMA 1 =GAMMA 2 I 770 

FIAM1=FLAM2 I 780 

48 ETA1=ETA2 l 7g 0 

DELZN=DELZN2/( TEMP 1-T SAVE ) /UHET 1800 

50 DZN=DZN+DELZN 1810 

TWR ( J ) =TDMY 1820 

DELZNWf J) =DELZN 1820 

FLWRITt J)=FLAM2 

55 CONTINUE 1 04 O 

IF(JAT)57,57,320 10 5O 

320 DLAHT= 1 • O+AHT I860 

DELT8= (DZN*DELT6*DLAHT )**BETAD 18 70 

GO TO 324 1880 

57 DELT8=(FLTN*DZN*DELT6/2.0)**BETAD 1890 

324 DELTA=DELTA2*DELT8 1900 

60 DIAT=DI+FLTN*(DELTA+DLCC) 

67 ZN=DZN/DIAT 2030 


IF(LAOTR-l) 446,446,448 
446 IF(GSRLQ) 450,450,452 

450 ZBG 1 = 0.7903*ZN*PINLT*PINLT*PTUBT*DI**4 

ZB I G= SQRT(ZBG1/(FMAS*RGAS*TINLT*VISC*FF) ) 
GO TO 69 

452 ZBG 1 = 0.7903*ZN*RHAIN*PDRPT*DI**4 
ZB I G= SORT (ZBG1/FMAS/VISC ) 

GO TO 69 
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448 ZBG1=FMAS23*ZN 

Z BG2=ZBG1**(9.0/14.0> 

68 ZB I G = Z BG2*ZBG6*ZBG7 

69 FLATN=ZN/ZBIG 

DEBUG LAOTR » ZBG1 . ZB I G , Z N , PDRPT »DI 
NB I G=FLATN 

B1GN=FLDAT (NBIG1+1.0 
IFlFLATN-BlGN+l.O) 71,71,73 
71 B I GN = B I GN-1 « 0 
73 VNE=VNE1/(BIGN*DI2) 

I F ( LADTR-1 ) 456,456,454 
456 H I = 4.36*TCGAS/DI 
GO TO 458 

454 IF(GLMR)212,212,210 
210 HI2=(HI1*VNE*DI )**0.4 
GO TO 75 

212 HI2=(HI 1*VNE*DI )»*0.8 
75 H I = H I 5 *H I 2 /D I 
458 D I AT 36 = 3600, 0/D I AT 

DEBUG HI ,HI2 ,VNE *BIGN 

BRANCH AND COMPUTE UO 


2050 
<2 060 
2070 

2080 

2090 

2100 

2110 

2120 


2140 

2150 

2160 

2170 


IF(QLMR) 520,520,503 
503 DAOD I = DIAT/DI 

U001 = DA0DI/(3600,0*HI ) 

UOLN= ALOGIDAODI ) 

U002= DI AT*UOLN/ ( THERKF+THERKF) 

UHET P= 1.0/(U001+U002) 

GO TO 522 

520 UHETP=DI AT36*HI*DI 
522 DIFFU=UHETP-UHFT 

T ESH= ( 'JHETP+UHFT ) *TESUH/2 • 0 

77 ADIFFU=ABS(DIFFU)-ABS(TESH) 

I F ( ADIFFUI65, 85,78 

78 KCNT=KCNT+1 

I F ( KCNT-20 ) 14,14,82 

82 WRITE (6,83) 

83 FORMAT ( 1H0,32H TROUBLE-LOOK AT STATEMENT NO. 78) 

GO TO 120 

85 FLEN=FLR*DI AT/2.0 
FLEN2=FLEN*FLEN 
FLDI=DI AT+2.0*FLEN 
TLITT2=TLITT1»FLEN2/FLAMRD 
TLITT=TLITT2*TZERC 
TLITTW=12.0*TLITT 
WF2 =WF 1*F LEN 
WT2 = ZN*WT 1 

I F(GSRLQ >445,445,435 
435 DCC=DI +DLCC+DLCC 
DCC2=DCC*DCC 
WCTl=PIE*ZN/4.0 
WCT2=RHAT*(DI AT*DIAT-DCC2 ) 

WCT3=RHAC*( DCC2— D I *D I ) 

WT=WCT1#(WCT 2+WCT 3 ) 

GO TO 87 

445 WT= PIE*ZN*(RHAT*DELTA*(DI+DLCC2 +DELTA ) +RHAC*DLCC* ( D I +DLCC ) ) 
DEBUG WT,DLCC, DELTA, ZN 

37 WF=WF2*ZN*TLITT 


2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2 340 
2350 
2360 

2380 

2390 

2400 

2410 

2420 

2430 

2440 


2460 
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on no r\ n 


PANEL=ZN*DIAT*( 1 .0+FLR) 2470 

REYN=VNE*DI*RHAIN/VISC 2480 

169 IF(GSRLQ) 180.180*170 2490 

170 P124=l. 0/24.0 2500 

P524=5. 0/24.0 2510 

P924=9. 0/24.0 2520 

DH1= PDRPH*RHAIN 2530 

DH2= BIGN*FLDI/DH1 2540 

DH4= DH2**P524 2550 

DH5=0. 174*VISC**P124 2560 

DH6=FMAS**P924 2570 

DH=DH4*DH5*DH6 2580 

WLQH2=BIGN*FLDI*RHAH 2590 

455 WCH1=2.0*PIE*BIGN*FLDI 2650 

WCH2=RHAH*DELTA* ( DH+DELT A+DLCC+DLCC ) 2660 

WCH3=RHAC*DLCC*( DH+DLCC ) 2670 

WLQH = WCH 1 * ( WCH2+WCH3 ) 2680 

173 WLQ1 2 = 2 • 0*BI GN*FLDI 2720 

WLQI3=WLOI2*DH*DH 2730 

WLQI4=WLQI3+DI2*ZN 2740 

WL0I=PIE*RHAIN*WLQI4/4.0 2750 

WL02 =WT + WF 2760 

174 WLO=WLQ2+WLOH+WLQI 2770 


DEBUG LAOTR*REYN»VNE,DI ,RHA I N ,V I SC ,B I GN , JL23 
I F ( JL23-1 ) 607,609,601 

609 I F ( REYN-2300 • 0 ) 648,610,610 
648 WRITE (6,646) 

GO TO 601 

610 WRITEt 6,612 ) REYN 
GO TO 120 

607 I F ( REYN-2300 • 0 ) 602,614,614 
602 JL23= 1 
LA0TR= 1 

WRITE (6,624) RFYN 
GO TO 605 
614 JL23 = 2 


601 IF(TZWR) 477,477,485 
477 DIAFLT=1.0/(DIAT*FLATN) 

DO 479 J= 1 , JEND 2800 

DELZWI J)=DELZNW( J ) #D I AFLT 2810 

479 CONTINUE 2820 

WRITE (6,467) (TWR(J) *J=1, JEND) 2830 

WRITE (6*469) (DELZW(J) ,J=1, JEND) 2840 

485 WRITE (6,175)WLQ*WLQH,WLQI ,WT ,WF,DH,DIAT .PANEL ,UHETP * DELTA, BIGN, 2850 

IZBIG.ZN.HI ,VNE,TLITTW,REYN 

175 FORMAT ( 1H ,5H WLQ= E 1 3 . 6 , 3X , 5HWLQH= E 1 3 . 6 , 3X , 5HWLQ I =E 1 3 .6 , 3X , 2870 

15H WT = E13.6,3X,5H WF = E 1 3 . 6 , 3X , 5H DH = E 13 . 6/ IX , 5HDI AT =E 1 3 . 6 , 3X , 2880 

26HPANEL=E13.6,2X,6HUHETP=E13.6,4X,6HDFLTA=E13.6,2X,5HBIGN=E13.6, 2890 

33X,5HZBIG=E13.6/1X,5H Z N= E 1 3 . 6 , 3 X , 5H H I =E 1 3 . 6 , 3X , 5H VNE=E13.6» 2900 

43X,6HTLITT=E13.6,3X,6H REYN=E13.6 ) 


I F ( JL23-2 ) 120,613,120 
613 JL2 3 = 3 

I F ( REYN-3000 . 0 ) 616,120,120 
616 LAOTR= 1 

WRITE (6,646) 
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n n 


GD TG 605 

180 FC1=PDRPH*RHAIN 

FC4=VISC/FMAS**0.6 

FC2 = FLD I /FC1 

WH1 =FMAS2/ D DR p H 

WH5=( 1.0/RHAIN)**0.2 

WH8*PIF*PHAH 

WH2 =WH 1 *F LD I 

AH3=1.405*FLOI 

AT1=PIE*DI AT 

FC3=FC2*BIGN 

FC5=FC3**0.2 

FC6=FC5*FC4 

FSM=0.0382*FC6**0.208 

WH3=WH2*FSM*BIGN 

WH4 =WH 3**0 • 2 

DNE3=WH3/RHAIN 

DNE=0.596*DNE3**0.2 

DTW1=RHAIN/RHAUT 

DTW2=DTWl**0.2 

DTW=DTW2*DNE 

WHH= BIGN*FLDI 

WIH= PIE*WHH*( (0.75*DNE +DLCC2 +DE LT A ) *DEL T A* RH AH +<0.75*DNE+ 
1 DLCC)*RHAC*0LCC) 

WGH = PIE*WHH*( <0.75*DTW +DLCC2 +DELTA ) *DELT A*RHAH 
1 (0.75*DTw+DLCC) *RHAC*DLCC) 

WH= WIH+WOH 
W= WH+WF+WT 

DEBUG WH,W*WOH*WHH,WlH*DELTA*BIGN 


2930 

2940 

2950 

2960 

2970 

3000 

3010 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3210 

3220 

3230 

3240 

3250 


AH1=WH4*WH5 
AH2=1 • 0+DTW2 
AH4=AH1*AH2*AH3 
95 AH= AH4*B I GN 
AT= AT1*ZN 
IF(JAT) 97*97,328 
328 KCNT = 0 

I F ( JAT6)330, 330*395 
330 AHT C= AH/ AT 

FAHT=AHT-AHTC 

I F ( AB5(FAHT)/FAHT ) 334*334,336 

334 NSG=-1 

GD TO 338 
336 NSG= 1 

338 IF(JATl) 340*340*341 

340 NSVl =NSG 
FAHTS1=FAHT 
AHT 51 = AHT 
JAT 1=1 

AHT1=0.25*AHT1 
AHT = AHT1 
GO TO 14 

341 IF( JAT2)342*342*360 

342 IF( NSG+NSV1 ) 344*343*344 

343 JAT 2=1 
FAHT S2 = FAHT 
AHT S2 = AHT 
GO TO 370 

344 I F ( JAT 4—6 ) 348,385*385 
348 JAT 4 = JAT 4+1 


3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

3460 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 
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n n on n n 


GO TO ( 350,352 ) ,JAT3 3560 

350 AHT1=C.25*AHT1 3570 

AHT = AHT 1 3580 

JAT 3 = 2 3590 

GO TO 359 3600 

352 AHT2 = AHT2 + 0 • 9 3610 

AHT = AHT2 3620 

JAT 3 = 1 3630 

359 UHETP=UHETRD 3690 

GO TO 380 3650 

360 IF(NSV1+NSG)362*365*362 3660 

362 NSV1=NSG 3670 

FAHT S1 = FAHT 3680 

AHT S1 = AHT 3690 

GO TO 370 3700 

365 FAHT S2 = FAHT 3710 

AHTS2=AHT 3720 

370 AHT SV= AHT 3730 

AHTDS= ( AHTS1-AHTS2 )/( FAHTS1-FAHTS2) 3790 

372 AHT=AHTS1-FAHTS1*AHTDS 3750 

ATST=0.5*(AHT+AHT5V)*FATST 3760 

I F ( ABS( AHT-AHTSV)-ATST >390,390,38 3770 

380 JAT 5 = JAT 5 + 1 3780 

I F ( JAT5-30) 19,19,382 3790 

382 WRITE (6,383 )JAT5,AHTC 3800 

GO TO 395 3810 

385 WRITE (6,387>JAT9,AHTC*NSV1 3820 

GO TO 395 3830 

390 JAT 6 = 1 3890 

GO TO 19 3850 

97 AHT = AH/ AT 3860 


DEBUG REYN »VNE »DI , RH A I N , V I SC , B I GN , JL2 3 , L AO TR 
395 I F ( JL2 3-1 ) 632,639,630 
639 I F ( REYN-2300.0) 650,636,636 
650 WRITE (6,696) 

GO TO 630 

636 WRITE(6,612) REYN 
GO TO 120 

632 I F ( REYN-2300, 0 ) 638,690,690 
638 JL2 3 = 1 
LAOTR = 1 

WRITE(6*629) REYN 
GO TO 605 
690 JL2 3= 2 


630 I F t TZWR)962,962,975 

962 DIAFLT=1.0/(DIAT*FLATN) 3880 

DO 965 J= 1 ♦ JEND 3890 

DELZWt J)=DELZNW(J)*DIAFLT 3900 

965 CONTINUE 3910 

WRITE (6,967) (TWR(J) ,J=1, JEND) 3920 

WRITE! 6,969) ( FLWR I T ( J ) , J=1 , JEND ) 

WRITE (6,969) (DELZW(J) ,J=1, JEND) 3930 

975 WRITE (6,100)AH,DNE*DTW,BIGN,WT,WF,WH,W»ZBIG»HI , AT , AHT ,Z N , DELTA , 3990 

ID I AT .PANEL ,UHETP , VNE , TL I TTW ,REYN 


I F ( JL23-2 ) 120,692,120 
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non no t* nn no 


642 JL23=3 

IF(REYN-3000.0) 644,120,120 
644 L AD TR= 1 

WRITE (6,646) 

GO TO 605 
120 CONTINUE 
526 CONTINUE 
525 CONTINUE 

530 I F ( READ-1 ,0 ) 8,2,215 


100 FORMAT ( 1H ,5H AH=E 1 3 . 6 , 3X » 5H DNE = E 1 3 , 6 ,3X , 5H DTW = E1 3 . 6 , 3X , 3960 

15HBIGN=E13.6,3X,5H WT=E1 3 . 6 , 3X , 5H WF=E13 . 6/ IX , 5H WH=E13.6,3X, 3970 

25H W=E13.6,3X,5HZBIG=E13.6,3X,5H H I =E13 .6 , 3X ,5H AT=E13.6,3X, 3980 

35H AHT=E13.6/1X,5H Z N=E 1 3 . 6 , 2X ,6HDEL T A=F 1 3 . 6 , 3X , 5HD I AT=E1 3 .6 , 3990 

42X,6HPANEL=E13.6,2X,6HUHETP=E13.6,3X,5H VNE = E 13 .6/ IX , 4000 

5 5HTLI T=E1 3.6,3X,5HREYN=E13.6 ) 

383 F 0RMAT ( 1H ,28HTR0UBLE - SEE STATE, NO. 382 ,4X , 5HJAT 5= I 3 , 4X , 4060 

15HAHTC=E12.5 ) 4070 

387 FORMAT ( 1H , 2 8HTR0UBLE - SEE STATE. NO. 344 ,4X , 5HJ AT4= I 3 , 4X , 4080 

15HAHTC=E12.5, 17) 4090 

9 FORMAT (11E7.5) 4100 

467 FORMAT ( 12F10. 2 ) 4110 

469 FORMAT ( 12F10. 4 ) 4120 

487 FDRMAT (//,1H0,35H GAS, LIQUID, LIQUID METAL RADIATOR/ 

1 , 1H0 »73H NOTE - THE COMPUTED TEMP . , L AMDA AND DEL Z PROFIL 4130 

1ES ARE PRINTED OUT BELOW/ 


2 1H0 »65H THE OUTPUT IS FOR TURBULENT FLOW EQUAS., UNLESS OTHERWISE 

3 STATED) 


612 FORMAT ( 1H , 8H I NST ABL E , 4X , 6HRE YN =E12.5) 

622 FORMAT ( 1H ,18HF0R LAMINAR EQUAS ., 2X ,6HREYN =E12.5) 

624 FORMAT! 1H ,20HFQR TURBULENT EQUAS . , 2 X , 6HRFYN =E12.5) 

646 F OR MAT (1H0.31H OUTPUT FOR LAMINAR FLOW EQUAS.) 

524 FORMAT ( 14, ( 10E7.5) ) 

662 FORMAT <1H0,/,2CH INPUT DATA - LIQUID) 

663 FORMAT (1H0,/,17H INPUT DATA - GAS) 

666 FDRMAT (1H0,/,26H IN°UT DATA - LIQUID METAL) 

END 4150 


IBFTC TABLE DECK 

SUBROUTINE TABLE 


DIMENSION TA8L( 13*36) 

COMMON TABL,FLR»ENC,ETATnT*SLPME,TUSF 


TUSE= 0.0 
6 FLSAVF=FLP 
13 IF(FNC)79,18,18 
18 I F ( FLR-10.0) 20,100,100 

20 IFIFLR-8. 0)310, 301, 301 
310 IF!FLR-6. 0)21 ,101,101 

21 I F ( FLP-4. 0)22, 102, 102 

22 IF1FLR-3. 0)311, 302, 302 


15 


311 IF(FLR-2. 0)23*103*103 
23 IF(FLR-1. 5)312*303, 303 

312 IF(FLR-1.0)71*104.104 
104 1 = 1 

DELL=0 • 5 
FLRT = 1 .0 
GO TO 60 
103 1 = 3 

DELL= 1*0 
FLRT=2.0 
GO TO 60 
102 1=5 

DELL = 2 * 0 
FLRT=4.0 
GO TO 60 
101 1=6 

DEL L = 2 • 0 
FLR T = 6 • 0 
GO TO 60 
303 1=2 

FLR T= 1 * 5 
DEL L = 0 • 5 
GO TO 60 
302 1=4 

FLR T = 3 • 0 
DELL= 1 • 0 
GO TO 60 
301 1=7 

FLRT=8.0 
DELL=2.0 
GO TO 60 

100 IF( FLR-15. )299, 299*300 
300 I F ( FLR-20. 0)99, 99*30 

30 IF(FLR-30.0)98,98,31 

31 IF{ FLR-50. 0)97*97*32 
299 1=8 

FLRT=10*0 
DELL = 5 • 0 
GO TO 60 
99 1=9 

DELL=5 • 0 
FLRT=15.0 
GO TO 60 
98 1=10 

DELL=10.0 
FLRT=20.0 
GO TO 60 
97 1=11 

DELL = 2 0*0 
FLRT = 30 »0 
GO TO 60 

32 1 = 12 

IF(FLR-1050. 0)40*41*41 
41 FLR = 10 50 
GO TO 40 
40 DELL=1000.0 
FLRT = 50*0 
GO TO 60 

60 IF( FNC-5. 0)67*70*200 
200 IF ( FNC— 20*0) 209*211 *206 
206 IF(FLR-50. 0)370*370*371 
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371 ET A= * 1 4 5 
GO TO 37 2 

370 ET A= *86 56778 5+FLR*(~. 1701 98 4+FLR* ( .0 1 2 71 25 3 7-FLP* . 000 1 91 947 ) ) 

372 WRITE (6*110) FNC 

HO FORMAT { 1H0*42HFNC OUT OF P ANGF-S EE SUBROUTINE TAELE- r NC=E12*5) 
GO TO 71 

209 FJ=21.+(FNC~5.0) 

J = F J 

F A J = J 

FRAC=FJ-FAJ 
K = 2 

GO TO 61 

67 FJ=F NC/O. 25+1.0 
J = F J 
FA J = J 

FR AC=F J-F A J 
K = 0 

61 DEL 1 =T ABL (I * J ) -T ABL t I +1 * J ) 

DEL2=T ABL ( I * J+l ) “TABL ( 1+1 ,J+ 1 ) 

FRAC1= (FLR-FLRT) /DELL 

ET A 1 = T ABL ( I ♦ J ) -DEL 1 *FR AC 1 
ET A2 = T ABL ( I ♦ J + l )-DEL2*FRACl 
ETA = ET A1-ETA2 
IF{ K>1 >62*72*210 

210 I F ( K— 2 )213*213*?14 

62 ETAT0T=ETA1-FRAC*ETA 
SLPME=-4*0*ETA 

flr=flsave 

GO TO 71 

70 J=20 
<=1 

GO TO 61 
72 ETAT0T=ETA2 

5LPME=”4.0*ETA 
FLR=FLSAVE 
GD TO 71 

213 ETATDT=ETA1-FRAC*ETA 
SLPME =— ET A 
FLR=FLSAVE 

GO TO 71 

214 ETATOT = ET A2 
SLPME = -ET A 
FLR = FL SAVE 
GO TO 71 

211 J = 3 5 
K = 3 

GO TO 61 
C 

c 

C TD SEE IF THE POINT ( FLR * FNC ) IS TABULATED IN THE TABLE 

C TUSE = 1* THEN POINT IS TABULATED 

C TUSE = 0 * THEN POINT IS NOT TABULATED 

C 

C 

C 

71 IF (FLR-1.0) 424*401*401 

401 IF ( FLR-5 0* 0 ) 402*402*424 

402 IF (FRAC1) 424*410,403 

403 IF (FR AC 1-1*0) 424*410*424 

410 IF ( FNC ) 424*411 *411 

411 IF (FNC-20.0) 412*412*424 


IT 



non on 


412 IF (FRAC) 413*423*413 

413 IF (FNC-5.0) 424,423*414 

414 I F ( FNC-20 • 0 ) 424,423*424 


423 TUSE= 1.0 
GO TO 425 

424 TUSE= 0.0 

425 RETURN 


79 ETA= 1 • 0 
RETURN 
END 
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SIBFTC GASLM L I ST ,RFF .DECK 

C SAULE - GAS. LIQUID. LIQUID METAL RADIATOR OOIO 

DIMENSION TWR< ICO) .DELZNWf 100 ) .DELZW ( 1 00 ) . BCDMY ( 1 2 ) .FXEQI 1000 ) » 

1 FXAB(IOO) 

COMMON FXEQ.FXAB.FLR.FXGAB.JI NTG . MESH , TST . DTST , TS . T EMP2 . FL AMTY ♦ 

1 FLAMC.ETAA1 , GAM A 1 . ET AA2 .C-AMA2 , FLAM , FWR I TE . FWRTSV 

READ IS. 490) (BCDMY(J) , J= 1 » 12 ) 

READ (5.134) FNMESH . FMESH .TESUH , FATST , FJ I NTG , TST , DTST ,FWR I TE 
WRITEI6.490) (BCDMY(J) ,J = 1.12) 


WRITE (6,487) 0120 

215 READ ( 5.3 JSTEPH. ALPHA, BETA, EMACH .RHAP.VAVG.PNAN, FCCLU.TDAY, 0130 

1E d SL»ATICK»GSRLQ «QLMR ,TZWR 0140 

NMESH=RNMESH 0150 


2 PIE= 3.1415926 
FWR T SV= FWR I TE 
JINTC= FJINTG 
MESH= FMESH 

0200 


BRANCH - GAS OR L I QU I D , L I QU I D METAL RADIATOR 0210 

0220 

IF( GSRLQ) 160, 160,150 0230 

150 READ ( 5,9)SCVEL»FMAS,VISC,UHETRD.CPH,TINLT,TEXIT,RHAIN,RHAF. 0240 

1 RHAH,RHAT,DDRDT,PDRRH,THERKF,TCGAS,TS.RHAC 0250 

IF (GLMR) 664,664,665 

664 WRITE (6,662) 

GO TO 667 

665 WRITE (6,666) 

667 WRITE <6,4)SCVEL»FMAS,VISC,UHETRD,CPH,TINLT»TEXIT, RHAIN.RHAF, 0260 

1RHAH.RHAT .PDRPT »PDRPH,THERKF,TCGAS,TS ,RHAC ,STEPH» ALPHA, BETA, 0270 

2E MACH, RHAP, V AVG, P NAN.FCCLU.TD AY, EPSL, ATI CK, FNMESH, FMESH, 

3 T ESUH.FAT ST, FJINTG, TST, DTST, FWRITE, GSR LQ,QLMR,TZWR 

419 TCLG=TCGAS 0350 

RHAL=RHAIN 0360 

155 FF=1.C 0370 

GO TO 165 0380 

160 READ ( 5,9)SCVEL»FMAS,VISC»UHETRD,CPH,TINLT,TEXIT, RGAS,PINLT» 0390 

1 RHAF.RHAH.RHAT »PTUBT,PTUBH,THERKF,TCGAS,TS,RHAC 
WRITE (6,663) 

WRITE ( 6 ,4 ) SCVEL »FMAS ,VISC»UHETRD,CPH ,T INLT ,TEX IT , RGAS » p I NLT , 0410 

1RHAF,RHAH,RHAT,dTUBT,dTUBH,THERKF , TCGAS, T S.RH AC, STEPH, ALPHA, BETA, 

1 EYACH, RH AP, V AVG, D NAN, RCCLU, TD AY, EPSL, AT I CK, FNMESH, FMFSH, TESUH, 

3 FATST.FJINTG, TST, DTST, FWR I TE , GSRLQ ,QLMR , JZWR 

PDR°T=°I NLT*PTUBT 0450 

PDRPH= D I NLT f-PTUBH 0460 

PEX IT=PI NLT- (PDRPT+2.0*PDRPH) 0470 

RHA I N = P I NLT/ < RGAS*T I NLT ) 0480 

163 RHAUT=DFXIT/(RGAS*TEXIT) 0490 

TINTD=1.0/TINLT 0500 

TEXTD=1.0/TEXIT 0510 

TINTD2=TINTD**2 0520 

TEXTD2=TEXTD**2 0530 

FF2=TEXTD2-TINTD2 0540 

FF3=TEXTD2*TEXTD-TI NTD2*TINTD 0550 
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FF=1.5*FF2/(FF3*TINLT) 

c 

165 FLTN= 2.0 

C 

DELT = ( T I NLT -TEX I T ) /ENMESH 
TEMPRD=TINLT-DELT/2.0 
H 1 1 = RHA I N /V I SC 
H13=VISC*CPH/TCGAS 
IF(QLMR)204,204,201 
201 H I 4 = H I 3**0 • 4 

HI5=0.625*TCGAS*HI4 
GO TO 206 
204 H I 4 = H I 3**0 • 3 

HI5=0.023*TCGAS*HI4 
206 DELZN1=FMAS*CPH*DELT 

DELZN2=3600.0*DELZN1/PIE 
FMAS2=FVAS*FMAS 
DELT1= (62.45*RHAP/RHAT)**0.5 
DELT2= ( VAVG/SCVEL )** ( 2. 0/3.0) 

DELT 3= ( 0.6747E-04/RHAP )** (1 .0/3 .0 ) 

BETAD=1.0/(3.0*BETA) 

DELT4=-TDAY*ALPHA/AL0G(PNAN) 

DELT 5=2.0/ (3. 0*EMACH*BETA+2 .0 ) 

DELT6=PIE*DELT4*DELT5 

DELTA1=DELT1*DFLT2*DELT3 

5 DELTA2=2.0*ATICK*FCCLU*DELTA1 
FMAS23=20.3/FMAS 
ZBG3=RHAIN*PDRPT/FF 
ZBG4=Z8G3**( 5.0/14.0) 

ZBG 5 = ( 1.0/VI SC) **(1.0/14.0) 

6 ZBG7=ZBG4*ZBG5 
VNE1=4.0*FMAS/(PIE*RHAIN) 

WT1 =P I E*RHAT 

TS4 = TS**4 

7 WF1 =2. 0*RHAF 

8 READ (5*9) D I , FLR » FL AMRD , RE AD 
WRITE (6.101DI »FLR.FLAMRD 

10 F ORMAT ( 1H0 » 5H DI=E13.6,3X» 

15H FLR=E13.6*3X*7HFLAMRD=E13.6) 

C 

C 

IF( RHAC )460. 460.461 

460 DLCC = 0 . 0 
GO TO 464 

461 D I 04 = . 04*D I 

IFIDI04-. 00125) 417,415.415 
415 DLCC = D I 04 
GO TO 464 
417 DLCC=. 00125 
C 
C 

464 DLCC2= DLCC+DLCC 
JL23= 0 
LAOTR= 2 

605 IF ( GSRLQ) 306.306,310 
306 JAT=1 
JAT 1 =0 
JAT 2 =0 
JAT 3 = 2 
JAT 4 = 1 
JAT 5 = 0 


0560 


0580 

0590 

0600 

0610 

0620 

0630 

0640 

0650 

0660 

0670 

0680 

0690 

0700 

0710 

0720 

0730 

0740 

0750 

0760 

0770 

0780 

0790 

0800 

0810 

0820 

0830 

0840 

0850 

0860 

0870 


0910 


0310 

0320 

0340 


0950 

0960 

0970 

0980 

0990 
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non 


* 

J AT 6 = 0 

1 000 


AHT = 0*3 

1010 


AHT1=C.3 

1020 

cc 

o 

r*' 

AH T 2 = 0 • 3 

1030 


GO TO 312 

1 040 

310 

JAT = 0 

1 050 

312 

SE D = 2.0*STEPH*EPSL/PIE 

CALL FXX ( FXEQ* C XAB,FXGA8*FLR*MESH, J I NTG * C WR I T E ) 
TL I TT 1 = 2.0*STE D H*EPSL/THER<F 



D I 2 = D I *D I 

1110 

11 

ZBG6=DI**(12. 0/7.0) 



UHETP=UHFTPD 

1150 


K.CN T = 0 

1160 

14 

UHET=UHETP 

SEPLR= 5EP*( 1.0+FLP) /UHET 
DZN = 0.0 
FLAMC= FLAMRD 
GAMA1= 5EPLP 
ET A A1 = 1. 

FLAMT Y= 0.0 
TEMP2= TEMPRD 
CALL TOJ(TDMY) 

1170 

25 

TSAVE=TDMY 

1390 


TZERC=TDMY**3 

DELZN=DELZN2/ (TEMPRD-TSAVE) /UHET 

1400 


DZN=DELZN 

1420 


TWR(1)=TDMY 

1430 


DELZNW ( 1 ) =DELZN 

1440 


JEND=NVESH 
FLAMC= 0. 

1450 

30 

DO 55 J = 2 * JEND 
GAMA1= GAMA2 
ETAA1= ETAA2 
TEMP2= TEMP2-DELT 
FLA^TY= FLAW/TDMY**3 
CALL TOJ(TOMY) 

DELZN=DELZN2/ (TEMP2-TDMY ) /UHET 

1460 

50 

dzn=dzn+delzn 

1810 


TWR ( J ) =TDMY 

1820 


DELZNW ( J ) =DELZN 

1830 

55 

CONTINUE 

1840 


I F ( JAT >57,57,320 

1850 

320 

DLAHT = 1 • O+AHT 

1860 


DELT8= ( DZN*DELT6*DLAHT)**BETAD 

1870 


GO TO 324 

1880 

57 

DELT8= (DZN*DELT6 >**3ETAD 


324 

DELTA=DELTA2*DELT8 

1900 

60 

DIAT=DI+2.0*(DELTA+DLCC) 


67 

ZN=DZN/DIAT 
DEBUG DLCC*DI 

I F ( LAQTR-1 ) 446*446*448 

2030 

446 

I F ( GSRLO) 450*450*452 


450 

Z BG1 = 0.7903*ZN*PINLT*PI NLT*PTUBT*DI**4 
ZB I G= SORT ( ZBG1 / ( FMAS*RGAS*T I NLT*VI SC*FF ) ) 
GO TO 69 


452 

Z BG1 = 0.7903*ZN*RHAIN*PDRPT*DI**4 
ZB I G= S0PT(ZBG1/FMAS/VISC) 

GO TO 69 
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448 

ZBG1=FMAS23*ZN 



Z BG2 = 2 BG1 ** (9.0/14.0) 

2050 

68 

ZBIG=ZBG2*ZBG6*ZBG7 

2060 

69 

flatn=zn/zbig 

2070 


nbig=elatn 

2080 


BIGN=FL0AT(NBIG)+1.0 

2090 


IF(FLATN-BIGN+1.0) 71,71,73 

2100 

71 

B I GN = B I GN-1 . 0 

2110 

73 

VNE=VNE1/ (BIGN*DI2) 

2120 



IF(LAOTR-l) 456,456,454 


456 

H I = 4 . 36*TCGAS/D I 



GO TO 458 


454 

I F ( QLMR >212,212,210 


210 

HI2=(HI1*VNE*DI ) **0 . 4 

2140 


GO TO 75 

2150 

212 

HI2=(HI1*VNE*DI )**0.8 

2160 

75 

HI=HI5*HI2/DI 

2170 

458 

D I AT 36 = 3600. O/DI AT 



branch and compute uo 


IF(OLMR) 520 * 520 * 503 
503 DAOD I = D I AT/D I 

U001 = DAOD I / ( 3600 . 0*H I ) 

UOLN= ALOG(DAODI) 

U002 = DIAT*UOLN/ t THERKF + THERKF) 
UHFT D = 1.0/(U001+U002) 

GO TO 522 


520 UHETP=DIAT36*HI*DI 2190 

522 DIFFU=UHETP-UHFT 2200 

TESH=(UHETP+UHFT)*TESUH/2.0 2210 

77 ADIFFIJ = ABS(D JFFU)-ABS(TESH) 2220 

IF(ADIFFU)85,85,78 2230 

78 KCNT=KCNT + 1 2240 

I F ( KCNT-20 ) 14,14,82 2250 

82 WRITE (6,83) 2260 

83 FORMAT(1HO,32H TROUBLE-LOOK AT STATEMENT NO. 78) 2270 

GO TO 120 2280 

85 FLEN=FLR*DIAT/2.0 2290 

FLEN2=FLEN*FLEN 2300 

FLDI=DIAT+2.0*FLEN 2310 

TLITT2=TLITT1*FLEN2/FLAMRD 2320 

TLITT=TLITT2*TZERC 2330 

TLITTW=12.0*TLITT 2340 

WF2 =WF 1*F LEN 2350 

WT2 =ZN*WT 1 2360 

IF(GSRLO) 445,445,435 

435 DCC=DI +DLCC+DLCC 2380 

DCC2=DCC*DCC 2390 

WCTl=PIE*ZN/4.0 2400 

WCT2=RHAT*(DIAT*DIAT-DCC2 ) 2410 

WCT3=RHAC* ( DCC2-D I *D I ) 2420 

WT=WCT 1* ( WCT2+WCT3 ) 2430 

GO TO 87 2440 

445 WT= PIE*ZN»(RHAT*DELTA*(DI+DLCC2 +DELTA ) +RHAC*DLCC* ( D I +DLCC ) ) 

87 WF=WF2*ZN*TL ITT 2460 

PANEL=ZN*DI AT* ( 1 ,0+FLR ) 2470 
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REYN=VNE*DI*RHAIN/VISC 2480 

169 IF( GSRLO) 180. 180.170 2490 

170 P124=l .0/24.0 2500 

P524=5. 0/24.0 2510 

P924=9. 0/24.0 2520 

DH1 = PDRPH*RHAIN 2530 

DH2 = BIGN*FLDI /DH1 2540 

DH4 = DH2** D 5 24 2556 

DH5=0 . 1 74*V I SC**P 124 2560 

DH6=FMAS**P924 2570 

DH=DH4*DH5*DH6 2580 

WL0H2=BIGN*FLDI*RHAH 2590 

455 WCH1=2.0*PIE*BIGN*F|_DI 2650 

WCH2=RHAH*DELTA* ( DH+DELT A+DLCC+DLCC ) 2660 

WCH3=RHAC*DLCC*( DH+DLCC ) 2670 

WLGH=WCH1 * ( WCH2+WCH3 ) 2680 

173 WLQI2=2.0*BIGN*FLDI 2720 

WL0I3=WL012*DH*DH 2730 

WLQI4=WLQI3+DI2*ZN 2740 

WLQI=PIE*RHAI N*WLQI 4/4.0 2750 

WLQ2 =WT + WF 2760 

174 WL0=WL02+WL0H+WL0I 2770 

DEBUG WT.DLCC.DFLTA.ZN 


DEBUG LAOTR.REYN.VNE.DI . RHA I N , V I SC . B I GN . JL 23 
I F ( JL23— 1 ) 607,609,601 

609 IFtREYN-2300.0) 648,610,610 
648 WRITE (6,646) 

GO TO 601 

610 WRI TE ( 6 ,612 ) REYN 
GO TO 120 

607 IF ( REYN-2300.0) 602,614,614 
602 JL23= 1 
LAOTR= 1 

WRITE (6,624) REYN 
GO TO 605 
614 JL2 3 = 2 


601 IF(TZWR) 477,477,485 
477 DIAFLT=1.0/(DIAT*FLATN) 

DO 479 J= 1 , JEND 2800 

DELZWf J)=DELZNW(J)*DIAFLT 2810 

479 CONTINUE 2820 

WRITE (6,467) (TWR(J) ,J = 1, JEND) 2830 

WRITE (6,469) (DELZW(J) ,J=1, JEND) 2840 

485 WRITE (6,175)WLQ,WLQH,WLQI,WT,WF,DH,DIAT,PANEL»UHETP, DELTA, BIGN, 2850 

IZBIG.ZN.HI ,VNE,TLITTW,REYN 

175 FORMAT ( lH ,5H WLQ=E 1 3 . 6 , 3X , 5HWLQH=E13 . 6 , 3X , 5HWLQ I =E1 3 .6 , 3X , 2870 

15H WT=E13.6,3X,5H WF=E 1 3 . 6 , 3X , 5H DH=E13 .6/ 1 X » 5HDI AT=E1 3 .6 ,3X » 2880 

26HPANEL = E13.6,2X ,6HUHETP = El 3 . 6 , 4X .6HDELT A = E1 3 .6 , 2X , 5HB IGN= E 1 3 .6 . 2890 

33X,5HZBIG = E13.6/1X,5H ZN = E1 3 .6 , 3 X , 5H H I =E1 3 .6 ,3X » 5H VNE = E13.6, 2900 

4 3X,6HTLITT=E13.6,3X,6H REYN=E13.6) 


I E ( JL23-2 ) 120,613,12' 

613 JL23= 3 

IF(REYN-3000.0) 616,120,120 
616 LAOTR= 1 

WRITE (6,646) 
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GO TO 605 

lflO FCl=DDR Dl -'*RHAIN 

FC4=V I SC /FVAS**^ • 6 

FC2=FLDI/FC1 

WH1=FMAS2/PDRPH 

WH5 = ( 1 ,0/RHAIN>**0.2 

WHS = P I E*RHAH 

WH2 =WH 1 *FLD I 

AH3=1.40S*PLDI 

ATl=PIF*OIAT 

FC3=FC2*BIGN 

FC5=FC3**0.2 

FC6=FC5*FC4 

FSM=0.0382*FC6**0.208 

WH3=WH2*FSM*BIGN 

WH4=WH3**0.2 

DNF3=WH3/RHAIN 

DNF=0.596*DNF3**0.2 

DTWl = RHAIN/RHA(JT 

DTW2=DTW1**0.2 

DTW=DTW2*DNE 

WHH= BIGN*FLDI 

W I H= PIE*WHH*( <0.75*DNE +DLCC2 +DELTA ) *DELTA*RHAH +<0.75*DNE+ 
1 DLCC)*RHAC*DI_CC) 

WOH=PIE*WHH*( <0.75*DTW +DLCC2 +DFLTA ) *DFLT A*RHAH 
1 (0.75*DTW+DLCC)*RHAC*DLCC) 

WH= WIH+WOH 


2 930 
2940 
2950 
2960 
2970 
3000 
3010 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3210 
3220 
3230 
3240 
3250 


W= WH+WF+WT 

AH 1 =WH4*WH5 3260 

AH2 = 1 . 0+DTW2 3270 

AH4=AH 1* AH 2* AH 3 3280 

95 AH=AH4*BIGN 3290 

AT=AT1*ZN 3300 

IF(JAT) 97,97,328 3310 

328 KCNT = 0 3320 

IF(JAT6)330,330,395 3330 

330 AHTC= AH/AT 3340 

FAHT = AHT-AHT C 3350 

I F( ABS ( FAHT) /FAHT ) 334,334,336 3360 

334 NSG=- 1 3370 

Gn Tn 338 3380 

336 NSG= 1 3390 

338 IF(JATl) 340,340*341 3400 

340 NSV1=NSG 3410 

FAHTS1 =FAHT 3420 

AHT S1 = AHT 3430 

JAT 1=1 3440 

AHT 1 = 0 . 2 5*AHT 1 3450 

AHT= AHT1 3460 

GO TO 14 3470 

341 IF(JAT2)342»342»360 3480 

342 IF(NSG+NSV1)344*343,344 3490 

343 JAT 2 = 1 3500 

F AHT S2 = F AHT 3510 

AHT S2 = AHT 3520 

GO TO 370 3530 

344 I F ( JAT4-6 ) 348,385,385 3540 

?48 JAT 4 = JAT 4 + 1 3550 

GO TO (350,352) ,JAT3 3560 

350 AHT1=0.25*AHT1 3570 

AHT = AHT1 3580 
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JAT 3 = 2 

3590 


GO TO 354 

3600 

352 

AHT2=AHT2+0.4 

3610 


AHT=AHT2 

3620 


JAT 3 = 1 

3630 

354 

UHETP=UHETRD 

3640 


GO TO 380 

3650 

360 

IF ( NSV1+NSG) 362 *365*362 

3660 

362 

NSV 1 = N SG 

3670 


FAHT 5 1 = FAHT 

3680 


AHT SI = AHT 

3690 


GO TO 370 

3700 

365 

FAHTS2=FAHT 

3710 


AHT S2 = AHT 

3720 

370 

AHT $V = AHT 

3730 


AHTDS=( AHTS1-AHT52 ) / ( F AHTS 1 -F AHTS 2 ) 

3740 

372 

AHT=AHTS1-FAHTS1*AHTD5 

3750 


AT5T=0.5* ( AHT+AHTSV)*FATST 

3760 


I F ( AB5( AHT-AHTSVJ-ATST ) 3 90*390*38 

3770 

380 

JAT 5 = JAT5 + 1 

3780 


I F { JAT5-30) 14*14.382 

3790 

382 

WRITE ( 6* 383 )JAT5 * AHTC 

3800 


GO TO 395 

3810 

385 

WRITE ( 6 » 387 ) JAT4 * AHTC *N 5V1 

3820 


GO TO 395 

3830 

390 

JAT 6=1 

3840 


GO TO 14 

3850 

97 

AHT = AH/ AT 

DEBUG WH*W,WOH,WHH,WlH,DELTA,BIGN 

3860 


DEBUG REYN *VNE *DI .RHAI N * VI SC. BIGN * JL23 *LA0TR 
395 I F t JL23-1 ) 632,634,630 
634 IF{ REYN-2300.0) 650,636,636 
650 WRITE (6,646) 

GO TO 630 

636 WR I TE ( 6 *61 2 ) REYN 
GO TO 120 

632 IF( REYN-2300.0) 638,640,640 
638 JL23= 1 
L A0TR= 1 

WR I TE ( 6 ,624 ) REYN 
GD TO 605 
640 JL23= 2 


630 IF(TZWR)462»462,475 
462 DIAFLT=1.0/(DIAT*FLATN) 

DO 465 J= 1 , JEND 
DELZW( J)=DELZNW< J)*DIAFLT 
465 CONTINUE 

WRITE (6,467) (TWR(J) ,J = 1, JEND) 

WRITE (6,469) (DELZW(J) ,J=1, JEND) 

475 WRITE (6»100)AH,DNE*DTW,BIGN,WT,WF,WH,W*ZBIG,HI , AT ,AHT ,Z N , DELTA , 
1DIAT,PANEL,UHETP,VNE*TLITTW,REYN 


3880 

3890 

3900 

3910 

3920 

3930 

3940 


I F ( JL23-2 ) 120,642,120 
642 JL23=3 
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IF( REYN-3000.0) 644,120*120 

644 LAO TR= 1 

WRITE (6,646) 

GO TO 605 
C 
C 

120 IF(READ-1.0)8,2,215 
134 FORMAT! 7E10. 3 ) 

3 FORMAT ( 6 El 2 . 5 ) 

4 FORMAT ( 8E1 5 . 5 ) 

490 FORMAT ( 12A6 ) 

487 FORMAT (//*1H0,35H GAS, LIQUID, LIQUID METAL RADIATOR/ 

1 , 1H0 *73H NOTE - THE COMPUTED TEMP . »LAMDA AND DEL Z PROFIL 
1ES ARE PRINTED OUT BELOW/ 

2 1H0 *65H THE OUTPUT IS FOR TURBULENT FLOW EQUAS., UNLESS OTHERWISE 

3 STATED) 

383 FORMAT ( 1H »28HTR0UBLE - SEE STATE. NO. 382 , 4X , 5H J AT5= I 3 , 4X , 

1 5HAHTC = E 12 . 5 ) 

387 FORMAT ( 1H , 28HTR0UBLE - SEE STATE. NO. 344 , 4X , 5H J AT4= I 3 , 4X , 
15HAHTC=E12 .5 , 17 ) 

9 FORMAT ( 1 1E7 . 5 ) 

467 FORMAT ( 12F10. 2 ) 

469 FORMAT ( 12F10. 4 ) 

100 FORMAT ( 1H ,5H AH=E1 3 . 6 , 3 X , 5H DNE = E 1 3 . 6 , 3X , 5H DTW = E 1 3 . 6 , 3X , 

15HBIGN=E13.6,3X,5H WT = E 1 3 . 6 , 3X , 5H WF = E13 .6/ IX , 5H WH=E13.6,3X, 
25H W=E13.6,3X*5HZBIG=E13.6*3X,5H H 1 =E13 .6 , 3X , 5H AT=E13.6,3X, 

35H AHT = E13.6/1X,5H Z N = E 1 3 . 6 , 2 X , 6HDEL T A = E 1 3 . 6 , 3 X , 5HD I AT = El 3 .6 , 
42X,6HPANEL = E13.6,2X,6HUHETP = E13.6,3X,5H VNE = E 13 . 6/ 1 X , 

5 5HTLI T=E13.6 ,3X , 5HREYN=E13 .6 ) 

612 FORMAT ( 1 H ,8 H I NSTABLE ,4X , 6HREYN =E12.5) 

622 FORMAT ( 1H ,18HF0R LAMINAR EQUAS ., 2X ,6HPEYN =F12.5) 

624 FORMAT! 1H ,2CHF0R TURBULENT EQUAS . , 2 X , 6HREYN =E12.5> 

646 FORMAT (1H0.31H OUTPUT FOR LAMINAR FLOW EQUAS.) 


662 

FORMAT 

X 

o 

<M 

+■ 

#> 

o 

X 

H 

INPUT 

DATA 

- LIQUID) 

663 

FORMAT 

( 1 HO * / * 1 7H 

INPUT 

DATA 

- GAS) 

666 

FORMAT 

END 

( 1H0*/*26H 

INPUT 

DATA 

- LIQUID METAL) 


4040 

0160 

0190 


4130 


4060 

4070 

4080 

4090 

4100 

4110 

4120 

3960 

3970 

3980 

3990, 

4000 ' 


4150 


$ I BFTC TO J LIST, REF, DECK 

SUBROUTINE TOJ(TDMY) 

FOR SAULE - GAS-L-LM 

DIMENSION FXEQt 1000 ) *FXAB( 100 ) ,FXA! 1000) *THA( 1000) , 

1 FXTHA! 1000) ,FTGA< 1000) 

COMMON FXEQ,FXAB*FLR ,FXGAB, J TNTG ,MESH , TST ,DTST ,TSA,TEMP2, 

1 FLAMTY »FLAMC * ETAA 1 *GAMA1 ,ETAA2 ,GAMA2 , 

2 FLAM,FWRITE,FWRTSV 
FWRTSV = FWRTSV-1 .0 

5 FJI NTG=JINTG 
MESH1= MESH+1 

7 FMESH=MESH 
DX= 1 • O/FMFSH 
FLR 1 = FLR+1.0 
KTDMY=0 

8 TDMY=TEMP2 


BEGIN TDMY I NTER AT I ON 
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c 

c 

c 

c 

c 


c 

c 

c 


TO 50LVE THE DIFFERENTIAL EQUA. USING SUB. 
AND SUB. FXX 


10 TDMY2=TDMY*TDMY 
THSA=TSA/TDMY 
FLAM=FLAMTY*TDMY**3+FLAMC 
12 THS A3=THSA**3 

THSA4=THSA3*THSA 
14 THSA5=1. 0— THSA4 
FLAMA5=FLAM*TH$A5 
16 FLAMA4=FLAM*THSA4 
DO 30 J=1,MESH1 

16 FXA ( J ) =— F L AMA5*FXE0 t J ) -FLAMA4 
30 CONTINUE 

I F ( FLAM ) 26.22,26 

22 DO 25 J= 1 »MESH1 
THA ( J ) =1 . 

25 CONTINUE 
GO TO 28 

26 CALL DEQ2(THA,FXA»FLAM»4. 0 .MESH, T ST » FWRT SV ) 

28 SLPA=(THA(2)-1.0)/DX 

TO COMPUTE ETA AND GAMMA 

I F ( FWRTSV >33.33,130 

130 WR I TE (6,132) DX , FLR 1 , TDMY »TDM Y2 »THSA , FLAM , THS A3 , THSA4 , 
1 THSA5 ,FLAMA5 ,FLAMA4 »SLPA 

132 FORMAT C 1H , 8H TO J- 133 , ( /8E15 .5 ) ) 

33 DO 35 J=1,MESH1 

34 FXTHA ( J ) = FXEQ(J)*THA(J) **4 

35 CONTINUE 

CALL FNTGRL t MESH1 »DX» FXTHA, FTGA ) 

ET1 A=2 . 0-THSA4 

37 ET2 A= ET1 A*FXGAB —FTGA { MESH1 ) 

IF(FLAM)38,39,38 

39 SLPMA=0.0 
ET A A2= 1. 

GAM AD= 0. 

GO TO 40 

38 SLPMA=SLP A/FLAM 

43 ET3 A= THSA5 +FLR* ( ET2A-SLPMA ) 

ETAA2= ET3A/THSA5/FLR1 

40 GAM A2= GAMA1*ET AA2/FTAA1 


C 

C 

c 

c 

c 

c 


END ETA AND GAMA CO M P. 


TO COMPUTE DERIVATIVE OF ETA WRT TDMY 


I F ( FLAM ) 44 ,62,44 
44 DRTHA1=4.0*THSA3*T SA/TDMY2 
TAFLR= THSA5*FLR1 
DXX=1.0/FJINTG 
DXX4= 4 . 0*DX X 
JDEL=MESH/JINTG 
JMESH=— JDEL/2+1 
I F ( FWR TSV)46,46,140 

140 WR I TE ( 6 ♦ 142 ) ET 1 A , ET2 A , E T3 A , FT AA2 .GAM A2 , 
1 DRTHA1,TAFLR,DXX,DXX4,SLPMA, 


DEQ2 
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1 ( FXEQ ( I ) * FXA ( I ) ,THA t I ) » FXTHA t I ) * 

1 THAI I ) ,FXTHA< I ) .1=98 ,101 ) 

142 FORMAT ( 1H , 8H TO J- 1 A 2 * ( /8 FI 5 . 5 ) ) 

46 ETAAD3=0. 


TO COMPUTE AN APPROXIMATION TO THE DERIVATIVE 
OF THETA WRT. TDMV 


DO 55 J=1*JINTG 
JME SH = JMF5H+JDFL 
TBARA = THA(JMESH) 

48 TBARA3 = TBARA**3 

DRTHA2 = DRTHA1*(DXX-FXAB(J) ) 

50 DRTHA= -DRTHA2/ (DXX4*TBARA3 ) 

ETAAD1=4.0*TBARA3* ( DXX-FXAB ( J ) ) 

52 ETAAD2= DRTHA*FT AAD1 

54 ETA AD3 = ET AAD3 + FT AAD2 
IF( FWRTSV155.55.150 

150 WRITEI6.152) JMESH , TB AR A , TB AR A3 . DR TH A2 . 

1 DRTHA.ETAAD1 . FTAAD2 . ETAAD3 
152 FORMAT ( 1H , 8H TOJ-152. 15. ( /8E15.5 ) ) 

55 CONTINUE 

ETAAD4= DRTHA1*FLR1*(1.0-ETAA2) 

58 ET A AD = ( ET A AD4 + FLR*ET A AD3 ) /TAFL.R 
60 GAM AD= GAMAl’FTAAD/^TAAl 
C 
C 

C END COMP. OF DERIVATIVE OF ETA WRT TDMY 

C 

C TO SOLVE FOR TDMY BY NEWTONS METHOD 

C 

62 TDM Y3= TDMY2*TDMY 
TDM Y4= TDMY3*TDMY 
TDMY34=4.0*TDMY3 
FNWT A= TDMY4-TSA**4 
64 FNWT = TDMY-TEMP2+GAMA2*FNWT A 

DFNWT A= GAMA2*TDMY34 + GAMAD*FNWT A 
66 DFNWT = 1 • O + DFNWT A 
TS AVE= TDMY 
TDMY= TDMY-FNWT /DFNWT 
68 TDTST= 0.5*ABS(TSAVE+TDMY)*DTST 
I F ( FWRTSV170. 70*160 

160 WRITEI6.162) KTDMY , ETAAD4 , ET AAD .GAMAD . 

1 TDM Y 3. TDM Y4 , TDMY 34, FNWTA, FNWT. DFNWT A, DFNWT , 

2 T SAVE, TDMY .TDTST 

162 FORMAT I 1H , 8H TO J- 1 60 , I 5 . ( / 8 El 5 . 5 ) ) 

70 I F ( ABS ( TSAVE-TDMY ) -TDTST 190.90.85 
85 KTDMY= KTDMY+1 

I F I KTDMY-30 1 10,10.87 

87 WRI TEI 6,88 1TDMY, TDTST .KTDMY 

88 FORMAT I 1H ,16H TROUBLE SUB TO J , 4X ♦ 5HTDMY = E 12 . 5 , 4X , E 1 2 . 5 , 4X , I 3 1 
90 RETURN 

END 


$ IBFTC FXX DECK 

SUBROUTINE FXX ( FXEQ , FXAB , FXGAB , FLR .MESH ,JI NTG 1 
C 

C FOR SAULE - GAS-L-LM AND SOLAR PROGRAMS 

DIMENSION FXEQI1000) »FXAB( 100) 
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3 FMESH=MESH 
FJINTG=JINTG 
5 DX= 1 « 0 /E^F 
MESH1=MESH+1 
X = — DX 

7 FLRD=1.0/FLR 

FLRDS= FLRD*FL p D 
FLRD2= FLRD+2.0 
9 FLRD5=0.5*FLRD 


TO OBTAIN FUNCT I ON ( X ) AT EACH MESH POINT 
FOR USE IN SUB DF02 

DO 30 J= 1 » MESH1 
X= X+DX 

16 FXA 1 =FLRD+X 

FX A2= FXA1*FXA1— FLRDS 
FXA= 1 .O-SORT (FXA2 ) /FXAl 
FXB 1 = FLRD2-X 

18 FXB2= S0RT(FXB1*FXB1-FLPDS) 

FXB= 1 .0-FXB2/FXB1 
20 FxEQ ( J ) = 0.5* ( FXA+FXB ) 

30 CONTINUE 
C 
C 

C TO OBTAIN JINT INTEGRALS OF 

C 

C INTEGRALS TO BE USED IN SUB 

C 

X = 0.0 
FX AB1 =0 . 0 
34 DX=1 .O/FJINTG 

FXGB1=FLRD2/FLRD 
FXGB2 = SQRT t FXGB1*FXGB1-1 .0) 

36 FXGB3=FXC-B2-ATAN(FXC-B2) 

DO 50 J=1,JINTG 
X=X+DX 

FXGA1= (FLRD+XI/FLRD 
FXGA2=SORT(FXGA1*FXGA1-1.0) 

38 FXGA3= FXGA2-ATANIFXGA2) 

FXGB4= ( FLPD2-X ) /FLRD 
FXGB5=SQRT(FXGB4*FXGB4-1 .0) 

40 FXGB6=FXGB5-ATAN ( FXGB5 I 
FXGB7=FXGB6-FXGB3 
FXAB2=-FLRD5* (FXGA3-FXGB7 ) 

FXABI J)=DX+FXAB2-FXAB1 
42 FXAB1=FXAB2 
50 CONTINUE 
C 

52 FXGAB=1.0+FXAB? 

55 RETURN 
END 


$ IBFTC DE02 DECK 

SUBROUTINE DE02 ( TH , FX , C .POW *MESH . TST ) 

C 

C FOR SAULE - GAS-L-LM AND SOLAR PROGRAMS 

DIMENSION TH< 1000 ) .FX ( 1000) . FD( 10001 »DG( 1000) »DGD ( 1000) 


FUNCTION (X) 
TOJ 


n n n non 


FMESH=MESH 
14 DX= 1 • 0 /FMESH 
DX2=DX*DX 
DGC=C*P0W*DX2 
16 D M1 =P0W-1 • 0 

FDC= ( —PM 1 )*C*DX2 
1R KCNT=0 

DO 20 J=1.MESH 
TH ( J ) = 1 .0 
20 CONTINUE 

22 DGD ( 1 ) =0. 

DG( 1 ) = 1 • 0 

23 FD( 1 ) = 1 . 0 

TO REDUCE MATRIX FROM TRI-DIAGONAL TO TWO DIAG 

DO 27 J=2.MESH 
THPOW = TH< J)**POW 

25 DG(J)= -DGD( J-l )-2.0-DGC*THP0W/TH( J) 

DGD ( J ) =1 . O/DG ( J ) 

FD1 =FDC*THP0W+DX2*FX ( J ) 

26 FD( J)=FD1-FD< J-l )/DG( J-l 1 

27 CONTINUE 

DGD < MESH )=1.0/(DG( MESH 1+1.0 ) 

28 TH ( ME SH+1 1=0.0 
MESH2=MESH+2 

30 KTST=0 

TO SOLVE FOR TH BY BACK SUBSTITUTION 

DO 35 JJ = 2 »MESH 
J=MESH2-JJ 
THSV = TH ( J 1 

31 TH( J)=DGD( J)*(FD( Jl-TH( J+l 1 1 

I F ( ABS ( THSV— TH ( J 1 ) -TST 1 3 5 ♦ 3 5 , 32 

32 KTST=KTST+1 

35 CONTINUE 

36 IF(KTST)43*43»37 

37 KCNT = KCNT + 1 
IF(KCNT-25)22»22»40 

40 WR I T E ( 6 * 4 1 1 KCNT >KTST 
43 TH ( MESH+1 1 =TH t MESH ) 

41 FORMAT ( 1H ,23H TROUBLE SEE SUBR. DEQ2 » ' X , 5HKC NT= I 3 . 4X , 5HKTST = 1 3 1 
45 RETURN 


COMPUTER PRINTOUT SHEET 
Description and Symbols 

This section contains a printout sheet from the electronic digital computer 
for an argon gas and an ether (ET-378) example. The sheet is typical for the 
inputs and outputs discussed in the analysis section of reference 1. 

The first group of numbers under the title Input Data represents the phys- 
ical and thermal properties of the working fluid, radiator material, and mete- 
oroid protection parameters. A brief explanation of these inputs, follows for 
a gas case as they are read from left to right. The symbols preceding these 
terms represent these inputs in the computational procedure of the computer 
program listing. 

First line : 

SCVEL sonic velocity, ft/sec 

TMAS total mass flow rate, lb/sec 

VTSC viscosity of gas, lb/(ft)(sec) 

UHETRD first guess of overall heat-transfer coefficient, Btu/(hr)(sq ft)(°R) 
CPH specific heat of gas, Btu/(lb)(°R) 

TINLT inlet gas temperature, °R 

TEXIT exit gas temperature, °R 

RGAS gas constant, ft-lb/(lb) (°R) 

Second line: 

PINLT inlet gas pressure, lb/sq ft 

RHAF fin density, lb/cu ft 

RHAH header density, lb/cu ft 

RHAT tube density, lb/cu ft 

FTUBT pressure drop fraction in tubes 

PTUBH pressure drop fraction in each header 

THERKF thermal conductivity of fin, Btu/(hr)(ft)(°R) 

TCGAS thermal conductivity of gas, Btu/(sec) (ft ) (°R) 
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Third line : 

TS equivalent sink temperature of space, °R 

RHAC liner density, Tb/cu ft 

STEEH Stefan-Boltzmann constant 

ALPHA constant in meteoroid mass distribution 

BETA constant in meteoroid mass distribution 

EMACH velocity ratio exponent in meteoroid protection equation 

RHAP meteoroid density, lb/cu ft 

VAVG average meteoroid velocity, ft /sec 

Fourth line: 

PNAN probability of no meteoroid penetration 

FCCLU orientation or occlusion factor 

TDAY radiator mission time, days 

EPSL emissivity of surface coating 

ATICK finite plate thickness and spalling correction factor 
FTEST error test in computation of surface temperature 

TEST error test in computation of surface temperature 

TESUH percent error in overall heat-transfer coefficient 
Fifth line: 

FNMESH number of elemental isothermal strips 

FATST percent error in header-area to tube-area ratio 

GSRLQ branch number (gas or liquid) 

QLMR branch number (liquid or liquid metal) 

TZWR control for printing of surface temperatures, conductance parameter, 
and length of elemental strip profiles 

The second group of input data is for a radiator that uses liquid ether 
(ET-378) as a working fluid. It is similar to the inputs of the gas example 
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except that instead of RGAS (gas constant), RHAIN (liquid density, lh/cu ft) is 
printed. The inlet gas pressure PINLT is not printed for the liquid example. 
Instead of FTUBT (pressure drop fraction in tubes) and PTUBH (pressure drop 
fraction in each header) FDRFF (pressure drop in tubes, lb/sq ft) and EDRFH 
(pressure drop in each header, lb/sq ft) are printed. The difference in inputs 
of branch numbers GSRLQ and QLMR should also be noted. For the gas example 
these inputs are 0. and 0.; for the liquid example, 0.10000E 01 and 0.; and for 
the liquid metal example (not printed), 0.10000E 01 and 0.10000E 01, respective- 
ly. 


These sets of inputs are for a simplified computational procedure when 
tabulated relations among effectiveness, conductance parameter, fin-tube pro- 
file ratio, and sink temperature are used. When a functional relation is used, 
additional inputs are included between inputs FATST and GSRLQ. 

FJIHTG used in approximation of integral in effectiveness computations 

TST absolute error for temperature ratio, 0^ 

O 

DTST percent error in surface temperature T 

FWRITE control for debug printing 

The group of numbers sifter the input data starts with parametric inputs. 

DI inside tube diameter, ft 

FLR fin-tube profile ratio 

FLAMED initial conductance parameter 

The next line represents the surface or base temperatures at the middle of 
each elemental strip (TWR(J)), followed by a line of corresponding conductance 
parameter (FLWRIT(j)) and length of each elemental strip (DELZW(J)). 

For the gas example, the following outputs are listed on the printout 


sheet , 

reading from left to right: 


AH 

total header area, sq ft 


DNE 

maximum inside diameter of 

inlet header, ft 

DTW 

maximum inside diameter of 

outlet header, ft 

BIGN 

number of tubes 


WT 

tube weight, lb 


WF 

fin weight, lb 


WH 

header weight, lb 
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W total radiator weighty lb 

ZBIG tube length, ft 

HI inside film heat-transfer coefficient, Btu/(sec)(sq ft)(°R) 

AT total tube area, sq ft 

AHT ratio of header area to tube area 

ZN total tube length, ft 

DELTA tube and header armor thickness, ft 

DIAT outside tube diameter, ft 

PANEL panel planform area, sq ft 

UHETP overall heat-transfer coefficient, Btu/(hr)(sq ft)(°R) 

VNE inlet fluid velocity, ft/sec 

TLIT fin thickness, in. 

REYN Reynolds number 

For the liquid example, most of the symbols are the same as for the gas 
example. The outputs and their symbols, that follow, are not used in the gas 
example . 

WLQ total radiator weight, lb 

WLQR header weight, lb 

WLQI liquid content weight, lb 

DH header inside diameter, ft 

TLITT fin thickness, in. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 8, 1965. 
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